Mon. Not. R. Astron. Soc. OOP. [TlfTTI C2009') Printed 1 April 2009 



(MN WI^ style file v2.2) 



A semi-automatic procedure for abundance determination 
of A- and F-type stars 

S. Hekker^*, Y. Fremat^, P. Lampens^, P. De Cat^, E. Niemczura^, 
O.L. Creevey^, J. Zorec^ 

^ Royal Observatory of Belgium, Ringlaan 3, 1180 Brussels, Belgium 

^Astronomical Institute, Wroclaw University, ul. Kopernika 11, 51-622 Wroclaw, Poland 
' ^ Instituto de Astrofisica de Canarias, C/ Via Lactea s/n, La Laguna 38205, Tenerife, Spain 

'^Institut d'Astrophysique de Paris, UMR7095 CNRS, Universite Pierre & Marie Curie, 98bis Bd. Arago, 75014 Paris, France 

, Accepted: Received; 



ABSTRACT 

A variety of physical processes leading to different types of pulsations and chemical 
compositions is observed among A- and F-type stars. To investigate the underlying 
mechanisms responsible for these processes in stars with similar locations in the H- 
R diagram, an accurate abundance determination is needed, among others. Here, we 
describe a semi-automatic procedure developed to determine chemical abundances 
of various elements ranging from helium to mercury for this type of stars. We test 
our procedure on synthetic spectra, demonstrating that our procedure provides abun- 
dances consistent with the input values, even when the stellar parameters are offset 
by reasonable observational errors. For a fast-rotating star such as Vega, our analy- 
sis is consistent with those carried out with other plane-parallel model atmospheres. 
Simulations show that the offsets from the input abundances increase for stars with 
low inclination angle of about 4°. For this inclination angle, we also show that the 
distribution of the iron abundance found in different regions is bimodal. Furthermore, 
the effect of rapid rotation can be seen in the peculiar behaviour of the H/3 line. 

Key words: stars: A/F-type stars - stars: fundamental parameters - stars: abun- 
dances - stars: individual: Vega 



1 INTRODUCTION 

A large variety of phenomena in terms of pulsations and 
chemical peculiarity is observed among main-sequence A 
and F stars. The underlying mechanisms of these different 
phenomena in stars with similar luminosity and temperature 
are not (yet) well understood. 

The p-mode pulsations in S Scuti stars are driven by the 
kappa (opacity) mechanism active in the second partial ion- 
isation zone of helium. For a long time it has been thought 
that these pul sations and c hemical peculiarity are mutually 
exclusive (e.g. lBregerlll970l ). The existence of chemical pecu- 
liarities as well as the suppression of pulsations due to a lack 
of helium in the partial ionisation zone can be ex plained by 
radiative diffusion (e.g. Michaud & Proffit3 [l993[ l. i.e. grav- 
itational settling of helium and levitation/settling of other 
elements in layers which are sufficiently stable against turbu- 
lent mixing. Notwithstanding this theoretical explanation, 
stars possessing both chemical peculiarities and pulsations 
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are observed. The co-existence of chemical peculiarity and 
pulsations can only be explained for more evolved stars, 
because the temperature range where He II partial ionisa- 
tion can occur sh ifts to deeper layers where some He is left 
(jCox et al.l I1979I ). In main-sequence stars the existence of 
both phenomena can not be fully understood in te rms of 
the radiative diffusion hypothesis. We refer to, e.g., iKurtj 
(2000) for an extensive overview. 

S Scuti and 7 Dor pulsators are found in similar re- 
gions of the H-R diagram, but the 7 Dor pulsation periods 
are an order of magnitude longer. This implies g-mode and 
not p-mode pulsations. The driving mechanism of the more 
recently discovered 7 Dor pulsators is convect ive blocking 
near t he base of their convective envelopes (e.g. Guzik et al.l 
|200gl; lLofflejl200i IWarner et all l2003l : iDupret et al] l2004l . 
200q). To improve the current models for 7 Dor stars, bet- 
ter constraints on the fundamental atmospheric parameters 
including chemical abundances are needed. A detailed abun- 
da nce analyses for abou t 20 7 Dor stars has been carried out 
bv lBruntt et all (120081 ). For this sample they find that the 
abundance pattern of the 7 Dor stars is not distinct from 
constant A- and F-type stars. Similar analyses for more stars 
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Table 1. Regions fitted to derive the projected rotational velocity 



wavelen 


jth range 
k 


4260 


- 4290 


4445 


- 4460 


4460 


- 4475 


4485 


- 4500 



Table 2. Spectral lines used to determine microturbulence for 
stars of various effective temperatures, 'all' indicates that a spec- 
tral line can be used for all considered temperature ranges. 



is needed to confirm this in a statistical sense. Only some 
stars show pul sation characteristics of both 7 Dor and 5 
Scuti type (e.g. lHenrv fc Feke]|l2005l ). Interestingly, most of 
these, so-called, hybrid stars appear to be metallic line (Am) 
stars. 

Furthermore, some very interesting spectroscopic mul- 
tiple syst ems with of A / F co mponents are known, such as 
DG L eo (jFremat et aLlbOOSal ) and 6"^ Tau l|Lampens et al.l 
I2OO7I ). These systems have at least one pulsating component 
and a non-pulsating or stable one. Further investigation of 
these systems in terms of chemical composition would al- 
low for a direct study of the relation between multiplicity, 
pulsations and chemical composition. 

Abundance determination is not only very interesting, 
as emphasised by the science cases described above, but 
also particularly challenging for A- and F-type stars, be- 
cause of the occurrence of rapid rotation and/or blended 
lines in these stars. So far, only a few chemical abun- 
dance analysis tools for these sta rs h ave been developed 



by, e . g. lErspamer fc NorthI (l2002bl) andlBruntt et all ||2002| . 



I2004l2008h which is based on IValenti fc Piskunovl (| 19961 ) 
With the spectra we have obtained for all science cases de- 
scribed above in mind, we independently developed a semi- 
automatic procedure to determine fundamental stellar pa- 
rameters (Section 2) and chemical abundances (Section 3), 
which we test on simulated data (Section 4). The main ad- 
vantage of such tests is that the input stellar parameters are 
well known a priory. These tests can provide an estimate of 
the error due to the procedure, and the uncertainties due 
to errors in the stellar parameters. Furthermore, gravita- 
tional darkening effects induced by fast rotation affect the 
abundance determination of stars seen nearly pole-on. We 
analysed Vega (Section 5) to compare results obtained with 
our method with literature values, and subsequently investi- 
gated the effects on abundance determi nation for simulated 
fast-rotating stars (|Fremat et al.|[2005bl l seen at several low 
inclination angles (Section 6). 



2 STELLAR PARAMETER DETERMINATION 

2.1 Model atmospheres and flux grids 

The radial velocity (RV), projected rotational velocity 
(vsmi), effective temperature (Tcff) and surface gravity 
(log g) are determined from spectra using a procedure 
partly based on the fitting of high-resolution spectroscopic 
data with synthetic spectra interpolated in a fiux grid. 
This grid was prepared by app lying the SYNSPEC pro- 
gramme (|Hubenv fc Lana 1 19951 . and references therein), 
and by adding opacities sources due to Rayleigh scatter- 
ing and to the H~ ions. All calculations were performed 



wavelength 

A 


element 


temperature range 
K 


4476.0 


Fe I 


< 8250 


4508.3 


Fe II 


all 


4515.5 


Fe II 


all 


4520.2 


Fe II 


all 


4558.7 


Cr II 


all 


4576.5 


Fe II 


< 8250 


4588.2 


Cr II 


> 7250 


4611.3 


Fe I 


< 7750 


4957.5 


Fe I 


< 9250 


5169.0 


Fe II 


> 7750 


5172.7 


Mg I 


all 



with the ATLAS9 LTE atmosphere models computed by 
ICasteUi fc Kurucd (|2003l l. The microturbulent velocity was 
fixed at 2 kms~'^, while a solar-type chemical composition 
(jCrevesse fc Sauvall 19981 ) was considered. 



2.2 Procedure 

The procedure we follow for stellar parameter determination 
(RV , ^'sin^, Teg, log g) c onsists of four consecutive steps (see 



also iFremat et all 2007h . First, RV is determined from the 
cross-correlation function (CCF) computed with a synthetic 
spectrum corresponding to the star's known spectral type. 
Then, v sin i is determined based on metallic line fitting in 
the sensitive regions indicated in Table [11 The fitting is per - 
formed with the Girfit computer code (jFremat et ahlliood ). 
where the Minuit Fortran package developed and applied by 
CERN is used to perfo rm least squ ares minimisation based 
on a simplex method (jPressI I2OO2I ) . All regions in Table [1] 
are fitted independently with the previously determined RV 
while log g, Tcff and vsini are free parameters. 

The stars we are concerned with are generally cooler 
than 8500 K, meaning that the hydrogen lines are sensi- 
tive to the effective temperature, and insensitive to Stark 
broadening and surface gravity variations. Therefore, we use 
Balmer lines to determine Tcff. For both A- and F-type stars, 
we use here the H/3 and Ha lines, since around these lines 
we can obtain reliable continuum normalisation. In early A- 
type stars, a reliable continuum could also be defined around 
the H7 and H(5 lines, but from tests on synthetic spectra we 
found that this is not the case for cooler stars, which are rel- 
atively more affected by metallic line blends. For this reason, 
and in order to conduct a homogeneous analysis of a sample 
of A and F stars, we do not take these lines into account. 
H/3 and Ha are fitted with the same packages as described 
above, with the previously obtained values for usini and 
RV, while log g is set to 4.0 in the first iteration. 

Finally, log g has been derived by combining the Hip- 
parcos parallax, apparent magnitude (reddening is not taken 
into account at the moment as we currently only have stars 
of 8*'' magnitude and brighter at our di sposal) and derived 
Tcff to obtain the stars' luminosity (jErspamer fc NorthI 
[2003). Tcff and luminosity values are then used to inter- 
polate the stellar mass and radius in the theoretical evo- 
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Figure 1. H/3 (top) and Ha (bottom) regions of tlie synthetic test spectra (black solid line) and the fitted (red dashed line) spectra for 
star I (left) and star II (right). (Colours only available in the online version.) 



Table 3. Stellar parameters of the two test spectra. The errors 
indicated for the computed values are the standard deviations of 
the values determined in different wavelength regions, systematic 
errors due to continuum normalisation have to be added. 





star I 


star I 


star II 


star II 




input 


computed 


input 


computed 


Teff [K] 


7500 


7367 ± 30 


7000 


6976 ± 14 


log g (c-g-s) 


3.5 




4.0 




vsini [kms~^] 


35.0 


35.8 ± 0.8 


32.1 


33.2 ± 0.4 


^micro [km S ^ ] 


3.70 


3.7 ± 0.5 


1.58 


2.2 ± 1.0 



lutionary tracks of ISchaller et al.l (|l992l ). which enables us 
to estimate the surface gravity. This value of log g is then 
adopted in a second iteration step to rederive the effective 
temperature. 

The results of applying the procedures described above 
to determine stellar parameters for synthetic test spectra 
and an observed spectrum of Vega are described in Section 
4, 6 and Section 5, respectively. 



3 ABUNDANCE ANALYSES 



to the wavelength region 4500 - 5500 A. These limits are set 
to avoid contamination with telluric lines, which occur at 
higher wavelength, and continuum normalisation problems. 
For wavelengths below 4500 A the high number of spectral 
lines in late A and cooler stars makes it hard to define the 
continuum for a large wavelength range, which might ham- 
per the abundance analysis. 

Model atmospheres and flux grids described in Section 
2.1 with stellar parameters determined with the procedure 
described in Section 2.2 are used for th e present analy- 
sis. Fu rthermore, the line list provided by iHubenv fc Lanz 
l|l995l ) is used with updated log gf values from the NISlLI 
database. 



Procedure 

1. Microturbulence parameter. Obtaining microturbu- 
lence {£,micro) is an important first step of an accurate abun- 
dance determination. W e adopt the method described by 
iGebran fc Monieij ([2003). They use the Mg II triplet at 4480 
A and a number of neighbouring unblended weak and mod- 
erately strong Fe II lines to fit for ^micro and rotational 
velocity while only allowing for small variations around the 
solar value. Here we use different spectral lines, listed in 



The abundance determination performed in the present 
work is based on a spectrum synthesis method and applied 
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Table (2] which are considered to be unblended and sensi- 
tive to ^micro in Certain temperature intervals. We fit these 
lines leaving ^micro as a free parameter and allow only for 
small changes in the abundances. The lines mentioned in Ta- 
ble [2] are all fitted separately, which is the approach we use 
throughout the whole procedure. To check the accuracy of 
the fit, and the sensitivity to a change in ^micro, the variance 
and a sensitivity parameter are calculated. This sensitivity 
parameter is the fractional change in fiux of each spectral 
line calculated for a synthetic fit with ^micro increased by 2 
kms~^ with respect to the best fit obtained. The 2 kms~^ 
level is chosen as an upper limit of the accuracy within which 
we want to obtain ^micro- The best fits are selected by eye 
and considering the variance and sensitivity parameter. The 
mean £,micro is computed and the standard deviation is used 
as an estimate of the error. 

2. Iron. We proceed by selecting isolated strong iron 
lines, i.e. lines with an equivalent width larger than 5 mA in 
a spectrum with the 'known' stellar parameters and so- 
lar abundance. Then, we search for the closest continuum 
around these lines in the spectrum and check whether there 
are no other strong lines in the wavelength interval selected. 
If this is not the case, the iron line is fitted, while ^micro can 
vary within its standard deviation around the mean value 
found in the previous step. For the sensitivity parameter of 
iron we take a 0.1 dex increase in the abundance compared 
to the best fit. The best fits are selected by eye, while also 
taking the variance and sensitivity into account. The aver- 
age value of the fitted ^micro is kept fixed throughout the 
remainder of the procedure, while the iron abundance is used 
as an input value and may vary within twice its standard 
deviation. 

3. Iron peak elements. Wc next fit the isolated iron 
peak elements, i.e. Ti, V, Cr, Mn, Co, Ni and Cu. We fit spec- 
tral regions where lines of these elements and possibly iron 
lines are strong. The variance and sensitivity parameter are 
calculated, and the best fitted spectral regions are selected. 
The mean values of the abundances found for each element 
are used as input for the abundance computation in spec- 
tral regions which contain only strong iron peak elements. 
From this analysis the abundance values for the iron peak 
elements are determined, which in turn are used to compute 
abundances for as many elements as possible ranging from 
helium to mercury. 

4. Final abundances. For the final abundance deter- 
mination, we consider all spectral lines with equivalent width 
^ 2 mA present in a spectrum with solar abundance, and 
a continuum is sought around these. All lines stronger than 
2 mA in the selected region are fitted, while the earlier de- 
termined ^micro is kept fixcd. The abundances of the iron 
peak elements obtained in the previous step are used (if 
present) as input, but these values can be changed within 
twice their standard deviation. For each element fitted in 
a certain region, the sensitivity parameter, as well as the 
variance, are computed. The worst fits are excluded and the 
mean abundances of all fitted elements are computed using 
the sensitivity parameter as weight. The standard deviation 
is used as an estimate of the error. We set the error to 0.1 
dex when the standard deviation is smaller than 0.1. For 
elements for which we only fit one spectral line the error is 
taken to be 0.2 dex. 

A last check is performed by inspecting the residuals 



Table 4. Input and com puted abunda nce differences with respect 
to the solar abundance IIGrevesse fc S auval 1998) for both star I 
and II. Nr indicates the number of spectral regions in which an 
element was fitted. 



el star I star I Nr star II star II Nr 



input computed input computed 



c 


-0.20 


0.05 ± 0.52 


15 


-0.11 


-0.01 ± 


0, 


.17 


17 


o 


-0.19 






0.03 


-0.02 ± 


0, 


.20 


1 


Na 


0.35 






-0.12 










Mg 


0.11 


0.08 ± 0.10 


5 


0.0 


-0.06 ± 


0, 


.10 


2 


Al 


0.28 






-0.21 










Si 


0.33 


0.20 ± 0.15 


9 


-0.09 


-0.16 ± 


0, 


.29 


13 


P 








-0.24 










S 


0.40 


0.15 ± 0.23 


2 




0.08 ± 


0, 


.10 


2 


Ca 


-0.11 


-0.16 ± 0.19 


7 


-0.07 


-0.07 ± 


0, 


.12 


7 


Sc 


-0.17 


0.00 ± 0.20 


1 


-0.01 


0.00 ± 


0, 


.10 


6 


Ti 


-0.09 


-0.10 ± 0.10 


31 


-0.13 


-0.16 ±0, 


.10 


33 


V 


U.DU 


n f^n -1-0 17 
u. ou m u. i / 


9 


u.uo 


0.04 ± 


0, 


1 A 


D 


Cr 


0.10 


0.11 ± 0.10 


37 


-0.04 


-0.05 ± 


0, 


.10 


57 


Mn 


-0.17 


-0.03 ± 0.10 


3 


-0.23 


-0.15 ± 


0, 


.10 


10 


Fe 


0.13 


0.04 ± 0.10 


93 


-0.10 


-0.18 ± 


0, 


.10 


114 


Co 


0.83 


0.54 ± 0.34 


3 


-0.01 


0.06 ± 


0, 


.10 


15 


Ni 


0.37 


0.34 ± 0.10 


36 


-0.16 


-0.15 ± 


0, 


.10 


43 


Cu 


0.32 


0.30 ± 0.20 


1 


-0.69 


-0.80 ± 


0, 


.13 


2 


Zn 


0.35 


0.53 ± 0.20 


1 


-0.35 










Sr 


0.82 






-0.07 










Y 


0.95 


0.92 ± 0.10 


4 


-0.01 


0.04 ± 


0, 


.27 


5 


Zr 


0.63 


0.75 ± 0.20 


1 


-0.27 










Ba 


1.18 


0.98 ± 0.27 


2 


0.39 










La 


0.94 


0.86 ± 0.20 


1 


0.20 


0.12 ± 


0, 


.20 


1 


Ce 


0.92 






0.11 


0.01 ± 


0, 


.10 


2 


Pr 








0.40 










Nd 


0.71 






-0.01 


0.08 ± 


0, 


.13 


4 


Sm 


0.83 






0.08 










Eu 


1.15 






-0.10 











between the observed and fitted spectra. In some cases, the 
residuals are larger for specific elements. The corresponding 
values are adjusted, to reduce the residuals and thus improve 
the final abundance determination. 



4 RESULTS FOR SYNTHETIC TEST 
SPECTRA 

To test the procedures described in the previous sections, we 
computed two synthetic test spectra with a signal to noise 
ratio of ~80 and a resolving power of ~40 000. For both spec- 
tra we determined the stellar parameters and abundances. 
The input and computed stellar parameters of the synthetic 
stars are given in Table |3] As described in Section 2.2, the 
surface gravities are computed from theoretical stellar evolu- 
tionary tracks, using the apparent magnitude and parallax. 
As the simulation and re-calculation would imply an exact 
inverse operation, we omitted such procedure here and only 
computed wsini and Tea, by taking into account the input 
log g value. We checked the influence of an offset in log g 
on the determination of Tcff by increasing / decreasing log 
g with 0.2. This induced changes in Tcff of less than twice 
the errors quoted in Table|3l Figure [T] shows the H/3 and Ha 
regions of the synthetic test spectra, together with the best 
fits we obtained for the stellar parameters. 
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C Mg Si S Ca So Ti V Cr Mn Fe Co Ni Cu Zn Y Zr Ba La 
elements 

Figure 2. Top: Region of input spectrum (black solid line) and fitted (red dots) spectrum for star I (left) and star II (right). Middle: 
Residuals between input and fitted spectra. Bottom: Input abundance (red) and computed abundances (black) of the individual elements 
determined using the method described in Section 3 (colours only in the online version). The thick error bars are the standard deviations, 
which serve as the error when they are larger than 0.1. The thin error bars are the inferred errors of 0.1 for an element fitted at several 
spectral ranges and 0.2 for an element fitted in only one spectral region. 



For the abundance analysis test, we used the input stel- 
lar parameters and the same atomic line list with which the 
spectra are computed. So, no errors due to stellar param- 
eters or inaccurate values in the atomic line list can occur 
(mainly errors in oscillator strength (log gf)). In this way 
we can really test the sensitivity of the procedure. In a sec- 
ond step (see below) the influence of offsets in the stellar 
parameters will be closely examined. 



For ^micro, we found 3.7 ± 0.5 kms~^ for star I and 
2.2 ± 1.0 kms~^ for star II, which is consistent with the 
input values provided in Table [S] Results of the abundance 
analyses are listed in Table [3] and shown in Fig. [2] together 
with the input and residual spectra. Throughout this work, 
the mentioned abundance values are al ways relative abun- 
dances with respect to the solar value l|Grevesse fc Sauvall 
il998iV For Ti, Cr, Fe and Ni most absorption lines with con- 
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C Mg Si S Ca Sc Ti V Cr Wn Fe Co Ni Cu Zn Y Zr Ba La Nd CO Wg Si S Ca Sc Ti V Cr Mn Fe Co Ni Cu Zn Y Zr Lo Ce Nd 




Figure 3. The computed abundances for different stellar parameters minus tlie input abundances (top) and minus the mean computed 
abundances (bottom) for star I (left) and star II (right). The black diamonds (solid black line) indicate the results of the computation 
with the correct stellar parameters of the synthetic test spectra. The coloured symbols are results from computations with different 
stellar parameters: red asterisks, green triangles, blue squares, light blue crosses for increased / decreased log g, Tcfi, Dsini and £,rn.icro 
respectively (colours only in the online version). The solid (dashed) lines connect values computed with increased (decreased) values of 
the stellar parameters with respect to the input values. The lines are only plotted for visual purposes. 



Table 5. £ micro for the tests with offsets in effective tempera- 
ture, surface gravity and rotational velocity. The input values for 
^micro are 3.7 and 1.58 kms~^ for test star I and II, respectively. 



star I star II 

Cmicro [km S ] £m,icro [km S ] 



Teff + 100 [K] 


4.2 


2.1 


TefF - 100 [K] 


3.8 


2.3 


log g + 0.1 (c.g.s) 


4.3 


1.3 


log g - 0.1 (c.g.s) 


4.2 


2.9 


V sin j + 5 [km s"-"-] 


4.3 


2.7 


V sini — 5 [km s^-'-] 


2.8 


0.2 



siderable strength are present and these elements are most 
important for the overall metallicity of the star. Apart from 
Mn and S for star I, abundances of all elements we com- 
puted are consistent, within their error-bars, with the input 
values. The mean difference in normalised flux between the 
input and fitted spectrum is -0.004 (input - fitted spectrum) 
for both stars with a standard deviation of 0.013 and 0.010 
for star I and II, respectively. 

These test spectra are also suitable to investigate the 



influence of an offset in effective temperature, surface gravity 
and rotational velocity. Therefore, we computed ^micro and 
abundances with a temperature of ± 100 K, leaving the 
surface gravity and rotational velocity at the correct value. 
Similar tests are performed for an increase / decrease of 0.1 
in surface gravity and 5 kms~^ in rotational velocity. 

The computed ^micro for the different tests are given in 
Table [5l The changes in S^micro due to changes in stellar pa- 
rameters are of the order of the error found for (^micro when 
we use the correct stellar parameters. Only for a decrease in 
wsini the value for S^micro decreases by about twice the er- 
ror. In order to test the influence of an offset in ^micro, the 
abundances are recomputed with (^micro ± 2 kms^^. The 
resulting abundances are shown in Fig. [3] 

When comparing the results of the tests described 
above, it becomes clear that the scatter of the computed 
abundance values around the input abundance of the syn- 
thetic test spectra is dominated by scatter inherent to the 
method (see top panels of Fig. [3} • The standard deviations 
of all computed abundances with respect to the input abun- 
dances are 0.17 dex and 0.12 dex for star I and II, respec- 
tively. The standard deviations around the mean of the com- 
puted abundances for the different stellar parameters are 
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Figure 4. Observed (black solid lino) and fitted (rod dots) speotra 
for Vega (top) with residuals (bottom) (colours only in the online 
version). 



0.085 dex and 0.10 dex for the respective stars (see bottom 
panels of Fig.|3]). From these tests we conclude that the accu- 
racy of the computed abundances is limited by the procedure 
itself and to a lesser extent by an inaccurate knowledge of 
the stellar parameters when uncertainties in these are equal 
to 100 K (Teff), 0.1 dex (log g), 5 kms"^ (vsini) and 2 
kms~^ {^micro) and when inaccuracies in the atomic line 
list can be neglected. Note that underestimating v sin i has 
the largest influence on ^micro and the abundances as seen 
from the light blue crosses connected with the dashed line 
in Fig. H 



5 RESULTS FOR VEGA 

As a subsequent test we analysed a spectrum of Vega (AO 
star with apparent magnitude of 0.03 mag in V) from the 
ELODIE database. Vega is often used as a reference star, al- 
though this is a fast-rotating star seen nearly pole-on, which 
influences the shape of the spectral lines. It is a bright star 
for which high signal-to-noise ratio spectra are available and 
for which several abundance analyses are already performed 
by other groups, providing comparison data for this test. 
The stellar parameters determined using our procedure are: 




Figure 5. Abundances of Vega determined using the method de- 
scribed in this paper (black). The coloured symbols are results 
from the literature: red dots, green triangles, blue diamonds an d 
light blue asterisks a r e resu lt s from lErspamer &: NorthI |[2002bl). 
Adcl man fc GulliveJ l ll990l) . lOiu et al.l 1I2OOII ') and lYoon et all 
II2OO8I ). respectively (colours only in the online version). 



Teff = 9560 ± 170 K, log g = 4.0 ± 0.3, vsini = 25 A ± 
2.3 kms"'^ and the radial velocity = -14.0 ± 1.0 kms~^, in 
goo d agreement with those det ermined from previous works 
(e.g. lErspamer fc Northll2002al ) . We obtained a ^micro of 2.4 
± 0.6 kms~^, and show a region of the best fitted spectrum, 
together with residuals in Fig. |4] The mean difference be- 
tween the observed and fitted spectrum of Vega is -0.0014 
(observed - fitted spectrum in relative intensity) with a stan- 
dard deviation of 0.007. The large feature at about 4860 A, 
H/3, could be due to deformations in the line caused by fast 
rotation, as explained in Section 6. 

Figure [5] compares the abundances we obtained 
for the different elements to those from _previous 



works, i.e., lAd ehnan fc Gulliver (1990), Qiu ct al. (2001!), 



lErspamer fcNorth (2QQ2b ) and Yoon et al. (2008). The last 
one is the only analysis taking gravitational darkening due 
to rotation into account, hence the difference in values of 
the stellar parameters at the poles and equator. The scat- 
ter in abundances between the different analyses is consid- 
erable. These differences might originate from the inclu- 
sion of the deformation due t o rotation or from the use 
of differe nt stellar para meters: 'A delman fc Gulliverl (|l990l ) 
and Ersp amer fc NorthI (2002b) us ed AT LAS model atmo- 
spheres available for Vega (|Kuruc3[l979l ) with T^ff = 9400 
K, log g = 3.9 and a fnncro of 2 and kms"'^, respec- 
tively, while iQiu et~aLl (|200ll ) used Teff = 9430 K, log g = 
3.95 and ^micro ~ 1.5 kms~^, which are somewhat cooler 
than our model atmosphere. Furthermore, for some ele- 
ments only a few lines can be fitted which make the fi- 
nal abun dance values very sensitive t o line parameters such 
as log gf. lErspamer fc NorthI l|2002bl ) argue that differences 
of some tenths of dex is not surprising for results from 
different authors, all using plane-parallel models. Our re- 
sults are in good agreeme nt with the abundances found by 
lErspamer fc NorthI ()2003h who used data obtained with the 
same instrument. Gravitional darkening due to fast rotation 
(such as in Vega) further introduces line profile deformations 
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Figure 6. Top; Fitted abundances for star III (left) and star IV (right). The thick error bars are the standard deviations, which serve 
as the error when they are larger than 0.1. The thin error bars arc the inferred errors of 0.1 for an element fitted at several spectral 
ranges and 0.2 for an element fitted in only one spectral region. Bottom: Residuals of the input spectra minus the fitted spectra for the 
full wavelength region. 



and biases on the chemical abundance determinations that 
arc more or less pronounced (|Yoon et al ] |2008l : lTaked~ et al.l 
[2008 ), depending on the line formation region (e.g., equator, 
poles or the whole surface). The offsets seen in Fig.[5]are con- 
sistent with their predictions. We further discuss the effects 
of fast rotation on chemical abundance determination in the 
next section. 



6 INFLUENCE OF FAST ROTATION ON 
PRESENT ANALYSES 

Here, we investigate what offsets in abundance we would ex- 
pect to find for fast-rotating stars seen nearly pole-on with 
our stellar parameter and abundance determination proce- 
dure. We simulated two spectra (star III and IV) with the 
same T^s and log g at the poles and the same rotational 
velocity, but with different incli nation angles. These simu- 
lations are based on model 4 of iTakeda et al.l (|2008l ). and 
have solar abundance. The input and computed parameters 
for both simulated stars are given in Table [B] 

The mean difference between the input and fitted spec- 
trum is 0.0003 and -0.001 for star III and star IV respec- 
tively with standard deviations of 0.0075 and 0.0045. The 
large spikes in the residuals are due to line deformations. 



The large feature in the residuals at about 4860 A is due 
to rotation and could serve as an indication of fast rotation 
when analysing real stars. The shape of the residual changes 
as a function of the difi'erence in stellar parameter values at 
the pole and at the equator. Nevertheless, the values we ob- 
tained for Tcff and log g are within the range of the input 
values. 

Results of the chemical abundance analyses are shown 
in Fig.|6] One notices that the abundances we derive are sen- 
sitive to the inclination angle, i.e., star III provides stronger 
deviations than star IV (both solar- like abundances). Grav- 
itational darkening effects due to rotation often lead to ap- 
parent chemical underabundances in AO stars. Star III in- 
deed shows underabundances ranging from to —0.25 for all 
considered elements, except for Al and S, which are found to 
be overabundant. It is interesting to note that we observe the 
same abundance pattern for Vega (Fig. [5]), although more 
pronounced. Also, the residuals in the H/3 line core are very 
similar for both Vega and star III. So, if Vega were at an 
inclination angle of about 4°, the gravitational darkening 
effects due to fast rotation could account for part of the 
chemical peculiarities of Vega. 

For star IV we recover the input abundances without 
any significant ofi^sets, which indicates that for stars with 
inclination angles > 7° we can recover the 'real' abundances. 
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Figure 7. Top: Fitted abundances for star IV with input abun- 
dances for Vega, indicated with red dots (colours only in the on- 
line version). The thick error bars are the standard deviations, 
which serve as the error when they are larger than 0.1. The thin 
error bars are the inferred errors of 0.1 for an element fitted at 
several spectral ranges and 0.2 for an element fitted in only one 
spectral region. Bottom: Residuals of the input spectra minus the 
fitted spectra for the full wavelength region. 



According to iTakeda et al.l l|2008l ) the inclination angle of 
Vega is ~ 7°, in which case we would be able to recover 
the 'real' abundances for this star. Computations of star IV, 
with the abundances we obtained for Vega confirm this (see 
Fig, 



m, 

lYoon et al.l (|2008l ) used a model for Vega with an in- 
clination angle of 4.54°, which is closer to our test star III. 
They mention that assigning a single abundance and (^micro 
for all lines of an element with low ionisation and excita- 
tion potentials, such as Fe I, is inadequate, because these 
elements are dominated by the contribution from the equa- 
torial region causing a double-horned shape in the spectral 
lines. A single abundance value for each element is however 
assumed in our procedure. To further investigate this, we 
study the distribution of the iron abundances fitted for dif- 
ferent spectral regions. The distributions for all four test 
stars are shown in Fig. [S] Clearly, star III with the lowest 
inclination angle has a bimodal distribution and thus the av- 
erage of the ab undance is incorrec t for most lines, confirming 
the findings of lYoon et al.l l|2008h . The bimodal effect is al- 



ready much less pronounced in star IV with an inclination 
angle of 7.4°. The distributions of the slowly rotating stars 
I and II are more symmetric and better centred. 



7 SUMMARY AND CONCLUSIONS 

In this work, we presented a semi-automatic procedure de- 
veloped to determine stellar parameters and abundances of 
elements from helium to mercury for A- and F-type stars. 
We performed various tests with this stand-alone procedure 
on synthetic spectra which had a resolving power of 40 000 
and a wavelength range of 4000-7000 A. 

Most abundance analyses are applied to a well-observed 
object like the Sun to verify and demonstrate its perfor- 
mance and dependability. In this study we investigate the 
accuracy of our analyses on synthetic test spectra. The main 
advantage of this is that we can perform the abundance anal- 
ysis with correct stellar parameters, i.e., the known Tefi, log 
g and vsini for which the test spectrum is computed and 
use similar log gf values for the "observation" as well as 
for the computations. This approach allows us to test the 
influence of erroneous stellar parameters on the abundance 
determination. From the analyses of these test spectra, we 
can conclude that we are able to obtain abundances which 
are consistent with the input values, within the errors. When 
"reasonable" changes (i.e. changes due to reasonable obser- 
vational errors) of the stellar parameters are considered, we 
show that the variance in the resulting abundances is smaller 
than the difference between the input and the computed val- 
ues. From this, we conclude that stellar parameters which 
have a reasonable offset from their 'real' values have only 
minor influence on the abundance determination. 

Our method is also tested on an observed spectrum 
of Vega, which is a fast rotator seen nearly pole-on, and 
compared with previous studies. The published results 
for Vega show that there are rather large differences be- 
tween the abundances from different studies. Our results 
are in good a greement with the abundances found by 
lErspamer fc North (2002U ) who used data obtained with 
the same instrument. 

Although we see only a small fraction of fast-rotating 
stars nearly pole-on, we might encounter one or more of 
them. We investigated this effect by simulating spectra 
which include fast rotation seen nearly pole-on. The com- 
puted abundances for the test star with an inclination angle 
of ~ 4° are offset from the input values, while this offset 
is negligible at an inclination angle of 7°. This is also re- 
flected in the bimodal distribution of the iron abundances 
obtained in different regions of the spectrum. This bimodal 
behaviour is likely due to the fact that lines with lower exci- 
tation potentials are more influenced by the rapidly rotating 
equatorial regions than lines with higher excitation poten- 
tials. For stars with an inclination angle of 7° the offsets 
from the input values are negligible and also the distribution 
of iron abundances obtained for different spectral regions is 
almost single-peaked. 

From our simulations we conclude that we should be 
able to recover the 'real' abundance values of Vega, if this 
star is indeed obser ved with an i n clinat ion angle of about 
7°, as suggested by iTakeda et al.l (|2008l ). Howe ver, if Vega 
has an inclination angle of about 4°, as used bv lYoon et al.l 
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Table 6. Stellar parameters of the two test spectra including fast-rotation. The errors indicated for the computed values are the standard 
deviations of the values determined in different wavelength regions. Q/Qc is the ratio of the rotational velocity to the break-up velocity. 
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Teff (pole) [K] 


9826 




9829 




log g (pole) (c.g.s) 


3.97 




3.97 




Teff (equator) [K] 


7971 




9408 




log g (equator) [cms^^] 


3.61 




3.89 




Teff [K] 




9133 ± 117 
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log g (c-g-s) 




3.80 ± 0.02 




3.95 ± 0.05 


V sini [km s"'^] 


15.7 


13.2 ± 2.4 
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13.6 ± 1.7 


inclination [deg] 


3.8 




7.39 
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0.44 
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2.9 ± 1.3 
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2.25 ± 1.2 
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Figure 8. Distributions of iron abundance obtained for the different regions of the spectrum for star I (left top), star II (right top), star 
III (left bottom) and star IV (right bottom). 



l|2008h . the influence of rotation should indeed be visible in 
the computed abundances. Note that the residual shape of 
H/3 of Vega is more similar to our simulations with the lowest 
inclination angle. 

The procedure for stellar parameter determination and 
abundance analysis, described in Sections 2 and 3, and 
tested on synthetic spectra and Vega, will now be applied to 
observations of A and F stars that we have at our disposal. 
This sample of stars contains both single stars and multiple 
systems, and these will allow us to study the relationship 



between pulsations, multiplicity and chemical composition. 
We have also obtained spectra of suspected 7 Dor stars, 
with the aim of both confirming their classification and ap- 
plying this abundance procedure. Results will be presented 
in subsequent publications. 
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